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Abstract 

We show that the ideal (nondissipative) form of the dynamical equations for the Lipps-Hemler 
formulation of the anelastic fluid model follow as Euler-Poincare equations, obtained from a con- 
strained Hamilton's principle expressed in the Eulerian fluid description. This establishes the math- 
ematical framework for the following properties of these anelastic equations: the Kelvin-Noether 
circulation theorem, conservation of potential vorticity on fluid parcels, and the Lie-Poisson Hamil- 
tonian formulation possessing conserved Casimirs, conserved domain integrated energy and an as- 
sociated variational principle satisfied by the equilibrium solutions. We then introduce a modified 
set of anelastic equations that represent the mean anelastic motion, averaged over subgrid scale 
rapid fluctuations, while preserving the mathematical properties of the Euler-Poincare framework. 

1 Introduction 

The Eulerian formulation of the action principle for an ideal fluid casts it into a form that is 
amenable to asymptotic expansions and thereby facilitates the creation and analysis of approximate 
fluid theories. Such an Eulerian action principle results whenever the general theory of reduction 
by symmetry groups is applied to Lagrangian systems, thereby yielding Euler-Poincare equa- 
tions, the Lagrangian analog of Lie-Poisson Hamiltonian equations, Marsden and Ratiu [1994]. 
This Euler-Poincare setting provides a shared mathematical structure for many problems in geo- 
physical fluid dynamics (GFD), with several benefits, both immediate (such as a systematic ap- 
proach to hierarchical modeling and versions of Kelvin's circulation theorem for these models) and 
longer term (e.g., structured multisymplectic integration algorithms). 

For example, by using the Euler-Poincare approach, Holm, Marsden and Ratiu [1998a,b], find 
that the action principles of a variety of incompressible fluid models for GFD are related by 
different levels of truncation of asymptotic expansions and velocity-pressure decompositions, as 
applied in Hamilton's principle for the unapproximated Euler equations of rotating stratified 
ideal incompressible fluid dynamics. This sequence of GFD models includes the Euler equations 
themselves, followed by their approximations, namely: Euler-Boussinesq equations (EB), primitive 
equations (PE), Hamiltonian balance equations (HBE), and generalized Lagrangian mean (GLM) 
equations. It also includes rotating shallow water equations (RSW), semigeostrophic equations 
(SG), and quasigeostrophic equations (QG). Thus, asymptotic expansions and velocity-pressure 
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decompositions of Hamilton's principle for the Euler equations describing the motion of a rotat- 
ing stratified ideal incompressible fluid are used in Holm, Marsden and Ratiu [1998a,b] to cast 
the standard models of GFD into Euler-Poincare form and thereby unify these descriptions and 
their properties at various levels of approximation. For related developments and additional struc- 
ture preserving approximations constructed from a similar viewpoint, see Allen and Holm [1996], 
Allen, Holm and Newbergcr [1998] and Holm and Zeitlin [1998]. 

Recently, Bannon [1996] reexamined the anelastic approximation for deep fluid convection and 
proposed an alternative form of the anelastic equations. This alternative model combines the 
results of Dutton and Fichtl [1969] and Lipps and Hemler [1982] to produce a hybrid theory that 
(1) conserves the domain integrated energy; (2) preserves potential vorticity on fluid parcels; and 
(3) accurately represents the acoustic adjustment process in Lamb's problem. The equations for 
a dry anelastic compressible fluid (atmosphere) rotating at angular frequency f2 under constant 
vertical gravitational acceleration gz take the following form Bannon [1996] 

du „ fp'\ gO' , . , 

^ +2f2xu = - v uJ + !: z < (L1) 

V-(p s u) = 0, (1.2) 
d(9 s + 8') 



0, (1.3) 

(1.4) 

(1.5) 



dt 

8' p' p' 

9 S p s gH P P. 



p' p' V 

Ps Ps T s ' 



In these equations, the fluid velocity is denoted u, the advective time derivative is d/dt = d/dt + 
u • V, and the state variables for this anelastic motion are: pressure p, density p, specific entropy 
8, and temperature T. These state variables consist of the sum of the base state (with subscript 
s) and a dynamic contribution denoted with a prime, as in 

d{x,y,z,t) = 8 s {z)+8'{x,y,z,t), (1.6) 

where 8' is the dynamic contribution to the specific entropy field. The base state is taken to satisfy 

d Ps DT n a dTTs ( , ^ 

~c[~ = ~9P s > Ps = p s RT s , C P 9s ~^ = ~ 9 > (I- 7 ) 

where ir s = T S /8 S . The constants R and C p are the ideal gas constant R and the specific heat at 
constant pressure C p for dry air. The density scale height is given by 1/H p = —p s (z)~ 1 dp s /dz. 
Given the base state functions satisfying relations (1.7), as well as the velocity u and the dynamic 
contributions p' and 8' at any time, the thermodynamic diagnostic relations (1.4) - (1.5) com- 
plete the description. The distinctions between this anelastic model and the traditional models 
Dutton and Fichtl [1969], Lipps and Hemler [1982] are discussed in detail by Bannon [1996]. For 
our purposes here, the important point is that the dynamical equations (1.1) - (1-3) agree in the 
formulations of both Lipps and Hemler [1982] and Bannon [1996] 

In this paper, we show that the dynamical equations (1.1) - (1-3) for the ideal (nondissipative) 
anelastic model follow as Euler-Poincare equations, obtained from a constrained Hamilton's princi- 
ple expressed in the Eulerian description. We then introduce a modified set of anelastic equations 
that represent the mean anelastic motion, averaged over subgrid scale rapid fluctuations, while 
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preserving the mathematical properties of the Euler-Poincare framework. Euler-Poincare equa- 
tions are the Lagrangian analog of Lie-Poisson Hamiltonian systems Marsden and Ratiu [1994], 
Holm, Marsden and Ratiu [1998a,b]. Among other things, the Euler-Poincare formulation of the 
anelastic fluid equations provides their Kelvin-Noether circulation theorem. This theorem is the 
basis for the conservation of anelastic potential vorticity on fluid parcels. Domain-integrated energy 
is also conserved and the relation of the Euler-Poincare equations to the Lie-Poisson Hamiltonian 
formulation of the anelastic dynamics is given by a Legendre transformation at the level of the 
vector fields satisfying the weighted divergenceless condition (1.2). The Casimir conservation laws 
for this Lie-Poisson Hamiltonian formulation provide a constrained-energy variational principle for 
the equilibrium solutions of anelastic dynamics and form a basis for determining their Lyapunov 
stability conditions, as done for the Euler-Boussinesq equations in Abarbanel et al. [1986]. (The 
Euler-Boussinesq equations form a special case of the anelastic model in which the base state is 
constant.) 

In related previous work, a two-dimensional study of the Hamiltonian structure of the Lipps- 
Hemler anelastic model was presented by Scinocca and Shepherd [1992], who also studied wave- 
activity conservation laws in the two-dimensional case. A canonical Hamiltonian formulation of 
the Lipps-Hemler anelastic model in three-dimensions was given in Bernardet [1995], Appendix 
A, for the Lagrangian fluid description of these equations. Perhaps not unexpectedly, the Lie- 
Poisson Hamiltonian formulation for the Eulerian fluid description of these equations agrees the 
canonical formulation of Bernardet [1995] and provides an alternative perspective. Indeed it must, 
as the Euler-Poincare theorem Holm, Marsden and Ratiu [1998a] proves for a class of Hamilton's 
principles that includes ideal continuum dynamics that the following four dynamical perspectives of 
fluid mechanics are equivalent: Hamilton's principle for the Lagrangian fluid description; the Euler- 
Lagrange equations in the Lagrangian fluid description; Hamilton's principle for the Eulerian fluid 
description with certain constrained variations, similar to those for reduced Lagrange d'Alembert 
equations; and the Euler-Poincare equations in the Eulerian fluid description. 

The methods of this paper are based on reduction of variational principles; that is, on La- 
grangian reduction. See Cendra et al. [1987, 1998a,b] and Marsden and Scheurle [1993a,b]), who 
also discuss systems with nonholonomic constraints. The latter has been demonstrated in the work 
of Bloch, Krishnaprasad, Marsden and Murray Bloch et al. [1996], who derived the reduced La- 
grange d'Alembert equations for such nonholonomic systems. Coupled with the methods of the 
present paper, these techniques for handling nonholonomic constraints would also be useful, if 
required, in continuum systems. In addition, it seems likely that applications of the techniques of 
multisymplectic geometry to numerical integrators associated with multisymplectic reduction will 
be exciting developments for the present setting; see Marsden et al. [1998] for a discussion of this 
approach. 

Organization of the Paper. In §2 we recall from Holm, Marsden and Ratiu [1998a] the results 
of the Euler-Poincare theorem for Lagrangians in continuum mechanics depending on advected 
parameters along with their associated Kelvin-Noether theorem and Lie-Poisson Hamiltonian for- 
mulation. These results establish the mathematical framework into which we place the dynamical 
equations for the anelastic model in §3. In §4 we introduce a modified set of anelastic equations that 
represent the mean anelastic motion, averaged over subgrid scale rapid fluctuations of amplitude 
a, while preserving the mathematical properties of the Euler-Poincare framework. 
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2 Applications of the Euler— Poincare Theorem in GFD 

Here we recall from Holm, Marsden and Ratiu [1998a] the statements of the Euler-Poincare equa- 
tions and their associated Kelvin-Noether theorem in the context of continuum mechanics and 
approximate models in geophysical fluid dynamics. 

The Euler-Poincare equations for a GFD Lagrangian L [u, D, b] involve fluid velocity u, buoy- 
ancy (or specific entropy) b and density D as functions of three dimensional space with coordinates 
x and time t. In vector notation, these equations are expressed as Holm, Marsden and Ratiu [1986, 
1998a,b,c] and Holm [1996], 

d 1 8L 1 5L„ ■ 1 , , 

dtD5u D 8vP D 5b 5D v ' 

or, equivalently, in "curl form" as, 

— — — - ux cur + V u H Vb = 0. (2.2) 

dt\D5uJ \D5uJ \ D5u 8DJ D Sb V ; 

The Euler-Poincare system is completed by including the auxiliary equations for advection of the 
buoyancy (or specific entropy) b, 

^ + u-V6 = 0, (2.3) 
and the continuity equation for the density D, 

^ + V-(Du) = 0. (2.4) 

For incompressible flows, one sets D = 1 in the continuity equation, so that V-u = 0. For anelastic 
flows, one sets D = p s (z) in the continuity equation with a prescribed stably stratified reference 
density profile p s (z), so that V • (p s (z)u) = 0. 

The Euler-Poincare motion equation in either form (2.1) or (2.2) results in the Kelvin- 
Noether circulation theorem, 

d f 1 5L 7 f 1 <5L 7 , 

dx.= - <t> — — Vb ■ dx. , (2.5) 



dt / 7t(u) D 5m J 7t(u) D 8b 

where the curve 7t(u) moves with the fluid velocity u. Then, by Stokes' theorem, the Euler- 
Poincare equations generate circulation of the quantity D~ 1 5L/5u whenever the gradients V6 and 
\7(D~ 1 5L/ 8b) are not collinear. 

Taking the curl of equation (2.2) and using advection of the buoyancy b and the continuity 
equation for the density D yields conservation of potential vorticity on fluid parcels, as 
expressed by 

^ + u ■ Vg = , where q = — V6 ■ curl ( — — J . (2.6) 
at D \D on J 

Consequently, the following domain integrated quantities are conserved, for any function $, 

C$ = J d 3 x D$(b,q), V$. (2.7) 
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The absence of explicit time dependence in the Lagrangian L [u, D,b] gives the conserved do- 
main integrated energy, via Noether's theorem for time translation invariance. This energy is 
easily calculated using the Legendre transform to be 

d 3 x (u-^J -L[u,D,b]. (2.8) 

When the Legendre transform is completed to express E [u,D,b] as H [m, D,b] with m = 5L/Su 
and 5H/5m = u, the Euler-Poincare system (2.1)-(2.4) may be expressed in Hamiltonian form 

^ = {u,H}, with ae[m,D,b], (2.9) 
and Lie-Poisson bracket given in Euclidean component form by 

{F,H}[m,D,b} (2.10) 



d 3 x 



5F 



5m 

6F, a . 5H 6F . . 5H \ 
h 6D^6m- + -Sb^Sm-j 



The conserved quantities C$ in (2.7) are then understood in the Lie-Poisson Hamiltonian for- 
mulation (2.9) - (2.10) of the Euler-Poincare system (2.1) - (2.4) as Casimirs that commute 
under the Lie-Poisson bracket (2.10) with any functional of [m, D, b]. The Casimirs also result via 
Noether's theorem from symmetry of the Hamilton's principle for the Euler-Poincare system un- 
der the "particle relabeling transformations" that leave invariant the Lagrangian L[u,D,b]. From 
the viewpoint of Noether's theorem, this particle relabeling symmetry corresponds to invariance 
of the Hamilton's principle for the Euler-Poincare equations under the transformation from the 
Lagrangian to the Eulerian fluid description, by pullback of the right action of the diffeomor- 
phism group on the configuration space of the Lagrangian fluid parcel positions and their veloci- 
ties. For full mathematical details, consult Marsden and Ratiu [1994]; Holm, Marsden and Ratiu 
[1998a,b,c]. 

The four properties (2.5)-(2.8) and the Lie-Poisson Hamiltonian formulation (2.9) - (2.10) 
of the Euler-Poincare equation (2.1) and its auxiliary equations (2.3) and (2.4) are desirable 
elements of approximate models for applications in geophysical fluid dynamics expressed in the 
variables [u,D,b]. Thus, the Euler-Poincare theory offers a unified framework in which to derive 
approximate GFD models that possess these properties: the Kelvin-Noether circulation theorem, 
conservation of potential vorticity on fluid parcels, and the Lie-Poisson Hamiltonian formulation 
with its associated conserved Casimirs and conserved domain integrated energy. Previous work of 
Holm, Marsden and Ratiu [1998a,b,c] has shown that many useful GFD approximations may be 
formulated as Euler-Poincare equations, whose shared properties thus follow from this underlying 
common framework. The aim of the next section of this paper is to cast the dynamical anelastic 
equations (1.1) - (1-3) into the Euler-Poincare framework, as well. 
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3 The dynamical anelastic equations are Euler— Poincare 
The Lagrangian 

In the Eulerian fluid representation, we consider Hamilton's principle for fluid motion in a three 
dimensional domain with action functional S = J dt L and Lagrangian L[u,D,b] given by 

L = J d 3 x D(^\u\ 2 + u-R(x)-gz-C p Ti s (z)b^ + p '(l-J?-j , (3.1) 

where D is the mass density and R(x), n s (z), and p s (z) are given functions of their arguments. 
This Lagrangian produces the following variations at fixed x and t 

ISL n/ . 5L „ . . 5L 1 D 

— — = u + R(x), — = -Cpir s (z) , — = 1 



D 5u 8b p v y ' 5p' p s (z) ' 

8L 1 v' 

— = -|u| 2 + u • R(x) -gz- C p n s (z)b - -t— . (3.2) 
oD 2 p s {z) 

Hence, from the Euclidean component formula (2.1) for Hamilton principles of this type, we find 
the motion equation for such a fluid in three dimensions, 

^-uxcurlR+ v(^\ + (g + C p b^jx = 0, (3.3) 

where curlR = 2f2(x) is the Coriolis parameter (i.e., twice the local angular rotation frequency). 
We use equation (1.7) to rewrite the last term in parentheses as 

dn s ( b\ 6'(x,t) , \ 

s+^^-^h^W' (3 - 4) 

in which we identify b as the total specific entropy, 

b = 9(x,t) = s (z) + 0'(x,t), (3.5) 

since each satisfies the scalar advection relation (1.3), cf. (2.3). Hence, from (3.3) and (3.4) we 
recover the anelastic motion equation (1.1), namely, 

!-„x2n(*)+vg)-£i = 0, (3.6) 

as the Euler-Poincare equation for the Lagrangian (3.1). Finally, we substitute the contraint D = 
p s (z) obtained from stationarity of the Lagrangian (3.1) with respect to variations in p' into the 
continuity equation (2.4) to find the anelastic divergence condition V- p s u = 0, i.e., equation (1.2). 
Preservation of this condition determines the dynamic pressure contribution, p', by solving the 
elliptic equation obtained by taking the divergence of the anelastic motion equation (3.8) after 
multiplying it by the base density p s (z), 

dp' 

— Ap' = g— h div (nonlinearity) . (3.7) 

The boundary condition for this elliptic equation is obtained from the normal component of the 
anelastic motion equation (3.8) evaluated on the boundary and using the boundary condition for 
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the velocity, e.g., that it has no normal component at the boundary, which yields a Neumann 
boundary condition for obtaining the pressure. 

Op 

— — — = gp (n • z) + n • (nonlinearity) . (3.8) 
On 

See the treatment in Bernardet [1995] for a discussion of alternative velocity boundary conditions 
for the Lipps-Hemler anelastic model. 

The Kelvin— Noether theorem 

From equation (2.5), the Kelvin-Noether circulation theorem corresponding to the anelastic motion 
equation (3.8) for an ideal anelastic fluid in three dimensions is, 

d_ 

dt 



(u + R)-dx = -rt C p 6Vir s (z) ■ dx , (3.9) 

7t(u) Jlt{u) 



where the curve 7t(u) moves with the anelastic fluid velocity u. By Stokes' theorem, the anelastic 
equations generate circulation of (u + R) around 7t(u) whenever the gradient of specific entropy 
9 has a horizontal component. Using advection of 9 and the anelastic divergence condition, one 
finds conservation of potential vorticity gAnei on fluid parcels, cf. equation (2.6), 

d^Anci + u _ y = q wn ere g Ane i = — rr • curl (u + R) . (3.10) 

dt p s {z) 

Consequently, the following domain integrated quantities are conserved, for any function $, cf. 
equation (2.7), 

C* = J d 3 x p s (z)$(9,q Ancl ), V$. (3.11) 

Energy conservation, Lie-Poisson Hamiltonian formulation and nonlin- 
ear Lyapunov stability analysis 

The conserved anelastic energy is easily calculated using the Legendre transform of the Lagrangian 
(3.1) to be 

E A nel = J d 3 X p s {z) Q|u| 2 + QZ + C p 7l s (z)9^j . (3.12) 

The corresponding Hamiltonian is (with b = 9) 

H Anel = J d 3 x (J^\m-DR\ 2 + Dgz + C p 7r s (z)Db\ + p' L^L - lj . (3.13) 

The Lie-Poisson bracket (2.10) now generates the dynamical anelastic equations (1.1) - (1-3) from 
this Hamiltonian according to equations (2.9). 

The canonical Hamiltonian formulation of the Lipps-Hemler dynamics due to Bernardet [1995], 
Appendix A, is based on the Hamiltonian in the Lagrangian fluid description, 

H canon = j d 3 a Q|±(a,t)| 2 + C p 7r s (^(a,t))&(a)^ . (3.14) 
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Transforming this Hamiltonian to the Eulerian fluid description yields, 

H A nei = Jd 3 x (J^\ m -DK\ 2 + C p 7r s (z)Dbj , (3.15) 

in which we again impose the anelastic density constraint explicitly, by using the pressure p' 
as a Lagrange multiplier. Applying the Lie-Poisson bracket (2.10) with this Hamiltonian yields 
the same dynamical anelastic equations as in (1.1) - (1.3), up to a redefinition of pressure to 
incorporate the gravitational acceleration. Thus, as guaranteed by the Euler-Poincare theorem 
Holm, Marsden and Ratiu [1998a] and the general theory of reduction Marsden and Ratiu [1994], 
the Lie-Poisson Hamiltonian formulation of the dynamical anelastic equations presented here in the 
Eulerian fluid description is equivalent to the canonical Hamiltonian formulation due to Bernardet 
[1995] in the Lagrangian fluid description. 

In the Eulerian fluid description we use the Casimir conserved quantities (2.7) to find the fol- 
lowing variational principle for anelastic equilibrium solutions: The equilibrium solutions 
of the dynamical anelastic equations occur at critical points of the sum H$, where 

H$ = H Ane i + C$ , (3.16) 

and 

C$ = J d 3 x D$(b,q), V$, where q = ^Vb ■ curl(m/£>). (3.17) 

Thus, the Casimir conservation laws for this Lie-Poisson Hamiltonian formulation of the three- 
dimensional anelastic equations in the Eulerian fluid description provide a constrained-energy 
variational principle for the equilibrium solutions of the dynamical anelastic equations and form 
a basis for determining their Lyapunov stability conditions, as done for the Euler-Boussinesq 
equations in Abarbanel et al. [1986]. The Euler-Boussinesq equations form a special case of the 
dynamical anelastic equations in which the base state is constant. Consequently, the analysis of 
the nonlinear Lyapunov stability conditions for the equilibrium solutions of the three-dimensional 
anelastic equations follows a similar procedure to that performed in Abarbanel et al. [1986] and 
produces a similar result, modulo the nonconstant base state. 

The pseudo-incompressible approximation (PIA) 

The PIA of Durran [1989] enhances the anelastic equations by allowing the influence of the base 
specific entropy field 9 s (z) on the mass balance. Additional studies of the PIA and compar- 
isons of its performance with the dynamics of the anelastic and Euler-Boussinesq models appear 
in Nance and Durran [1994], Lilly [1994]. 

As an Euler-Poincare system, the PIA equations arise from a modification of the anelastic 
Lagrangian in equation (4.1). Namely, we consider the Lagrangian for the PIA system given by 



-PIA 



d 3 x 



( i|u| 2 + u ■ R(x) -gz- C p n s (z) 6 



+p ( Ps (z)6 s (z)-De 



PIA constraint 
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This Lagrangian slightly modifies the volume constraint in the anelastic Lagrangian (4.1). The 
modified PIA constraint on D imposed by the pressure p requires 

D* = D6 = p s (z) 6 s (z). 

Hence, the fluid velocity u must satisfy a weighted incompressibility relation (pseudo-incompressibility) 

3D* 

— + V-(£>*u) = 0, V- ( Ps (z)6 s (z)u) = 0, (3.19) 

obtained via time independence of D*, and the advection equations for specific entropy 9 and 
density D, 

f)f) f)F) 

^ + u-W = 0, ^_ + V-( J Du)=0. (3.20) 

The PIA equations are the same as for the anelastic model except for the effects of the differently 
weighted incompressibility constraint in (3.19). 



4 The dynamical anelastic-alpha equations are Euler— Poincare 
The Lagrangian for the anelastic-alpha model 

Following Holm, Marsden and Ratiu [1998a,b], we introduce into the dynamical anelastic equations 
the effects of averaging over rapid fluctuations whose amplitude falls below the length scale denoted 
as a, by making the following modification of the Lagrangian in equation (3.1) for the anelastic 
model 



d x 



2 



D[ i|u| 2 + ^-|Vu| 2 + u ■ R(x) - gz - C p n s (z)b 

+ M D 



Ps(z) 



(4.1) 



where | Vu| 2 = & and we denote the other variables as before, in equation (3.1). This modified 
Lagrangian for the anelastic-alpha model produces the following variations at fixed x and t 

1 SL a 2 , ^ x _ . . SL ^ . . SL D 

u - — (Du k ) + R(x) , — = -C p 7i s (z) 



DSu D v ' v ' ' 5b " y 77 bp' p s {z 

|§ = i|u| 2 + ^|Vu| 2 + u • R(x) -gz- C p7 r s (z)b - ^- . 
oD 2 2 Ps\z) 



Hence, from the Euclidean component formula (2.1) for Hamilton principles of this type, the 
Euler-Poincare equation for the Lagrangian (4.1) is given by 

dv „ ■ „ „( v' 1, ,9 a 2 ,„ l0 \ q6' A , 

— + v j Vu J -u x curlR+ V — - -|u| 2 - — |Vu| 2 - ^—z = 0. (4.3) 



dt J \ps 2 1 1 2 J 9„ 

This is the motion equation for the anelastic-alpha fluid in three dimensions. One should compare 
this motion equation with equation (3.8) for the standard anelastic fluid and note that the advective 
time derivative is still defined as d/dt = d/dt + u- V, in terms of the fluid parcel transport velocity, 
u. We have used the following additional notation in equation (4.3), 

v = u-a 2 Au, and Au=—(Du k ) k , (4.4) 
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where the weighted Laplacian operator A is given by 

and subscript-comma index notation denotes partial spatial derivatives. The key observation about 
the anelastic-alpha model is that its transport velocity u is smoother relative to its momentum or 
circulation velocity v, by inversion of the Helmholtz operator, (1 — A), which depends upon the 
stratification profile of the base state. The curl form of equation (4.3) is 

dv , „ / p' 1 , l9 a 2 .„ ., \ q& . . . . 

— -u x curl(R + v) + V — — -|u| 2 — — - |Vu| 2 + u ■ v - ^-z = . (4.6) 
ot \p s 2 2 J 9 S 

As we shall see, the analysis of the anelastic-alpha model follows the same procedure as for the 
anaelastic model, with appropriate changes from u to v to account for the effects of the averaging 
over the subgrid (sub-a-length) scale rapid fluctuations. 



The Kelvin— Noether theorem for the anelastic-alpha model 

From equation (2.5), the Kelvin-Noether circulation theorem corresponding to the anelastic motion 
equation (3.8) for an ideal anelastic-alpha model in three dimensions is, 

v + R) -dx= - I C p 9Vir s {z) -dx, (4.7) 



where the curve 7f(u) follows the fluid parcels in moving with anelastic-alpha fluid velocity u. 
Thus, by Stokes' theorem, the anelastic-alpha equations generate circulation of (v + R) around 
7t(u) whenever the gradient of specific entropy 9 has a horizontal component. Using advection of 
9 and the anelastic divergence condition, one finds conservation of the anelastic-alpha potential 
vorticity qA-a on fluid parcels, which are transported with the anelastic velocity u, cf. equation 
(2.6), 

+ u ■ Vg A - a = , where q A ^ a = — ■ curl (v + R) . (4.8) 

Ot Ps(z) 

Consequently, the following domain integrated quantities are conserved, for any function $, cf. 
equation (2.7), 

C$ = J d 3 x p s (z)<S>(9,q A _ a ), V$. (4.9) 
Energy conservation for the anelastic-alpha model 

The conserved anelastic-alpha energy is calculated as before using the Legendre transform of the 
Lagrangian (4.1) and found to be 

Ea-o = j d 3 x Pa ( z )(l\ u \ 2 + ^\Vu\ 2 + gz+C p ir a (z)eY (4.10) 

Thus, the kinetic energy is augmented in the anelastic-alpha model by a term proportional to the 
squared amplitude of the velocity shear. Hence, the anelastic-alpha model costs energy to produce 
velocity shear and its solutions will tend to have smoother velocity profiles than those for the 
standard anelastic model. 
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The curl of equation (4.6) and the anelastic divergence condition V ■ p s u = yield an equation 
for anelastic-alpha vortex dynamics, 



The control on the L 2 norm ||Vu|| 2 afforded by the conserved energy in equation (4.10) should 
tend to moderate the vortex stretching term q ■ Vu in the vortex dynamics equation (4.11) and, 
thus, produce less violent, smoother and more coherent turbulent vortex dynamics than occurs for 
the standard anelastic equations. Numerical investigation of the anelastic-alpha system introduced 
here will be conducted and reported elsewhere. 
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